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ABSTRACT 


One candidate for a reusable upper stage to be carried by the Space 
Shuttle is an Aerosianeuvering Orbit-to-Orbit Shuttle (AMOOS). This con- 
cept uses the drag of the vehicle during a pass through the ataosphere 
rather than the propulsion system to slow the vehicle on a retuxu from 
a high energy orbit. In this paper, the nature and magnitude of the 
sensitivity of AMOOS to uncertainties in the properties of the atmosphere 
is shown. Various guidance schemes for correcting for the effects that 
the unpredictable variations in the atmosphere have on the trajectory are 
discussed. For the mission studied here, a payload retrieval from 
geosynchronous orbit with aerodynamic plane change, a linear feedback 
guidance scheme was developed. A relatively simple heuristic law was 
used to demonstrate the conce*>t. Using optimal control theory a feedback 
law was developed analytically. Testing with a large number of different 
atmospheres showed this law to be a feasible means of controlling the 
AM30S trajectory. Also, refinements to the technique offer promise of 
significant inq)rovement and these are discussed. 
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IHTRODUCnON 


One essential element of the Space Transportation being developed 
for future space operations is the reusable upper stage to be carried 
by the Space Shuttle. The current baseline vehicle is a purely pro- 
pulsive Orbit-to-Orbit Shuttle (OCS). An alternative candidate, the 
Aeromaneuvering Orbit-to-Orbit Shuttle (AMOOS), uses the drag of the 
vehicle during a pass through the atmosphere rather than the propulsion 
system to slow the vehicle on a return from high orbit. AMOOS offers 
much more payload capability than OOS and uses current rather than ad- 
vfiuiced technology. 

One of the key problem areas of AMOOS is its sensitivity to uncer- 
tainties in atmospheric properties and the developaent of a guidance law 
which will correct for these uncertainties and Insure that the vehicle 
flies the desired trajectory. This paper will describe the background 
of the AMOOS concept, compare it with the conventional OOS, and examine 
the guidance problem. The nature and the magnitude of the sensitivity 
of AMOOS to the state of the atmosphere will be discussed. This sensi- 
tivity W81S determined by linking a trajectory program with an atmosphere 
generator program which produces a large number of different atmospheres 
that statistically match actual atmosphere data. The data from a large 
number of trajectories provided clues for possible guidance schemes which 
use the aerodynamic capability of the vehicle. 

In partlculeur, a linear feedback guidance law will be developed and 
tested on a large nvunber of different atmospheres. First, a simple 
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heuristic lav vill he used to demonstrate the feedback technique. Then, 
using optimal control theory, a feedback lav vill be developed analytic- 
ally. Possible refinesients to the technique and areas of further study 
will also be proposed. 
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BACKGROUND 


To fully realize the potential of the Space Shuttle presently being 
developed under the direction of NASA, a reusable upper stage, or Orbit- 
to-Orbit Shuttle (OOS), becomes a necessary component of the total space 
transportation system. The OOS, placed in orbit by the Shuttle, would 
deliver and retrieve payloads from high energy orbits. Detailed studies 
of a purely propulsive OOS are presently underway using a baseline vehicle 
described in reference 1. The mission profile for such a vehicle and 
a schematic of its configuration are shown in figures 1 and 2. The maximum 
weight of the OOS is constrained by the maximum Shuttle payload capability 
of 29^8U kg. With this restriction, oven technology advances such as com- 
posite materials and a new improved propulsion system will not give OOS 
the capability to perform etll of the missions proposed in the payload 
mission model without the use of multiple Shuttle latanches. One study 
(r 2 f. 2) states that a significant number of missions will require two 
Shuttle/OOS launches, thus increasing the cost and complexity of those 
missions. This situation is fxirther aggravated by the high sensitivity 
of the performance of a purely propulsive OOS to changes in inert weight 
and engine specific impulse. One possible way to reduce this sensitivity 
and increase payload capability of a reusable upper stage may be to use 
the atmospheric forces rather than the engine thrust to brake the vehicle 
during the transfer from high orbit to Shuttle rendezvous altitude. This 
could eliminate about 1/k of the propulsive Av required for a round trip 
geosynchronous mission. 
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^ff&ibmmJTY OP THE 


The concept of using the atmosphere to alter a vehicle’s orhit has 
heen the subject of much study in the psist. Most of this work has been 
in the form of euialytical studies of synergetic plane change, which uses 
a combination of proptilsive and aerodynamic forces to change the inclina- 
tion of the vehicle’s orbital plane (refs. 3-12). The initied and final 
orbits of the vehicle are essentially the same except for inclination. 
These studies, limited to neai- earth orbit, indicate that the synergetic 
plane change can give substantial fuel savings for vehicles with hyper- 
sonic L/D more than 2. 

The technique of using the atmosphere to slow an orbital vehicle 
has been in use since the early days of space flight. However, aero- 
braking, as applied to full reentry from space requires that essentially 
all of vehicle’s energy be dissipated, while the application to a return 
to a low earth orbit frem synchronous orbit requires only a specific, 
precise amount of energy loss. 

In recent years, severad studies have investigated the potential of 
aerobraking for CX)S type vehicles. Multi-paiss trajectories were con- 
sidered to minimize aerodynamic heating. Guidance and heating constraints 
were determined (ref. 13) and an analytical technique for calculating 
sore trajectory parameters was developed (ref. lU) for midtlpass missions. 
The trajectory of a typical aerobraking mission is Illustrated in figure 
3. A later comprehensive study investigated the feasibility and practi- 
cality of the aerobraking mode for return from high orbit of em OOS with 
an aerobraking kit attached (ref. 15). Some of the conclusions of this 
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study vere that the aerobrahing node vas feasible, payload vas maximised 
by missions vith 25 to 33 aerobraUng passes, and that, contrary to con- 
ventional OOS, the aerobraking OOS could perfom the baseline round trip 
mission \ising current technology. 

The most recent activity in this area has been a feaislbllity study 
of a nev vehicle called the Aeromaneuvering Orbit-to-Orbit Shuttle (AMOOS) 
(refs. l6, 17). The results of this study include the identification of 
the most promising configurations for AMOOS, demonstration of potentially 
large payload gains, and identification of problem areas. One critical 
problem area vas the sensitivity of PiKXiS to uncertainties in the proper- 
ties, particularly density, of the atmosphere and the g\iidance of AMOOS 
through an unpredictable atmosphere. This problem and Its solution is 
the subject of this paper. 
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AMOOS 


The vehicle used in the present investigation vos derived in the 
AMOOS feasihility studies. This section provides a general description 
of AMOOS and the particular configuration studied herein using material 
froiD, references l6 and 17. 

The AMOOS mission vlll differ from that of the propulsive OOS only 
in the manner in which it achieves a phasing orbit with the Shuttle during 
the return phase. AIKX}S will be targeted to enter the earth's atmosphere 
for a prescribed number of passes. It vlll exit after the last pass 
with Just enovigh velocity to carry it to phasing orbit apogee. At apogee, 
a short burn is required to achieve phasing orbit perigee. Once this 
is done, the AMOOS mission is then identical to that of the propulsive 
OOS. 

The most promising configuration evolved as a compromise between 
high lift to maximize control authority, high drag to minimize TPS mass, 
and aerodynamic trim requirements. Also, volume and dimensioned constrednts 
imposed by both the Shuttle paylosid bay and potential payloads impacted 
the design. The configuration chosen is illustrated in figure U. Using 
em ablative thexmed-protection system on a one-pass mission yielded the 
lightest vehicle and thus the largest payload. The results of a cost 
euiedysir edso were favorable to the ablator (reference 17). 

The payload capability of AMDOS is shown in figure 5% compared to the 
propulsive OOS. The advantage over a purely propulsive OOS lies in missions 
to high orbits where the propellant savings as a result of aerobraking 
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are largest. Another important result of the feasibility studies is that 
AMOOS is indeed less sensitive than the OOS, to variations in specific 
impulse and structural mass. This is due to the smaller AV requirement 
which allovs for a higher inert mass. 

The principal disadvantages to AMOOS lie in the area of navigation 
and guidance. The combination of navigation errors and uncertai’^ties 
in the properties of the atmosphere could cause the vehicle either to 
reenter or fail tr dissipate enoiigh energy to effect Shuttle rendezvous. 
Furthermore, the time histories of most of the trajectory state para- 
meters are indistinguishable for a wide range of atmospheric density 
variations until the vehicle has very nearly reached its lowest point. 

This leaves little time to detect errors and make the necessary correc- 
tions (reference l8). 

Studies (references l6 and 17) have showr that the payload retricvsd 
mission, since it results in the heaviest vehicle at atmospheric entry, 
has the greatest impact on the design of vehicle systems required for 
atmospheric flight. Therefore, it will serve as the design mission for 
this paper. The mission is a return from geosynchronous orbit at 0® 
inclination to e, low altitude phasing orbit at 28-5° inclination. Most 
of the pl6me change is done propulsively during the first burn and the 
remainder is done using the AMOOS lift capability during the single pass 
throiigh the atmosphere. The mass of the vehicle at entry is 10338 kg and 
its aerodynamic coefficients are shown in Figure 6. 
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AMOOS GUIDANCE CONCEPTS 


This section is a summary of the potent** al guidance techniques for 
AMOOS, some of which depend on the +ype of mission. Hv..^ever, the basic 
concept is the same for all. That is, to correct for atmospheric uncer- 
tainties the drag of the vehicle is modvl^'ted by either raising or lowering 
the vehicle’s planned altitude profile. This can be done either with the 
main propulsion system or by controlling the vehicle aerodynamic ally. 

This study only considered the latter. A study of guidance concepts is 
contained in the AMOOS literature (reference 18) and some of the material 
in tnis document is repeated here. 

For the case when the original high energy orbit of AMOOS is in the 
same pleuie as that of the Shuttle, no plane change is required of the 
vehicle during either the orbital coast or the atmospheric portion of the 
return trajectory. In fact, while in the atmosphere, AMOOS must experi- 
ence no net lateral force to avoid a costly propulsive plane change 
maneuver after exiting ^‘le atmosphere. 

"’he most obvious wa** of controlling ailtitude when there is no lateral 
force requirement is direct modification of the lift vector via angle- 
of-attack. This technique can best be used if AMOOS flies inverted, that 
is, with a bank angle of l80°, to increase control authority. For in- 
st€ince, if the atmospheric density is lower than expected, AJ400S would 
increase the angle-of-attack to lower the altitude profile, thereby encoun- 
tering higher density. At the same time, the higher angle-of-attack results 
in a higher drag coefficient, contributing to a higher rate of decelerat.^on. 
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If the vehicle vere at 0^ hank angle* these effects iroold tend to cancel 
each other. Thus* using the Inverted node vDid.d decrease the sstount of 
angle-of^atteck through idJ,eh the vehicle must rotate to correct its path. 
A problem vitb the angle-of-attack control is that the moment of inertia 
of the vehicle in the longitudinal plane is quite hi^, making vehicle 
rotation in that plane difficult. 

The inertia is much less about the centerline, vhich leads to a con- 
sideration of baidc angle as the means of control. In this mode, the 
angle-of-attack remains constant and the vehicle would nominally fly 
at either 90® or 270® bank angle. Off nominal conditions could be cor- 
rected by appropriately modulating bank angle up or down depending on 
the nature of the conditions. The problem with this mode is that it 
would result in an orbit plane change. This could be eliminated by re- 
versinc, the bank angle. For instance, the vehicle would nominally be 
at a bank angle of 90 ® from the point of entry to perigee and 270® from 
perigee to exit. However, the vehicle would have to switch quickly, 
a difficult maneuver, to avoid excessive perturbation of the trajectory. 
Sank angle modulat ion was used in the development of a guidance law 
in reference 18. 

AMOOS could benefit from the advantages of the above modes by using 
some cond>ination of both. It could he that the trajectory could be 
divided into various segments during idiich one mode may be more applicable 
than the other. However, any scheme which would have the advantages, of 
both .jodes would also have the disadvantages and in addition would be 
more complex. 
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TUfi 
lOft 

Another method discussed In the AMOOS literature is a constant rota- 
tion of the vehicle about the velocity vector during the nontinal trajec- 
tory. This results in no net lateral force. In this case, off-nominal 
conditions are compensated for by varying the angular rate during appro- 
priate portions of the cycle in such a manner that the average rate remains 
unchanged over one revolution. The result is a net lift in a direction 
which may be controlled. The problems of this mode are that it is sen- 
sitive to the initieLL value of bank angle and it requires an almost con- 
stant firing of the vehicle's reaction control thrusters. 

Most of the missions to high energy orbits require a change in orbit 
inclination of 28.5°. On the return trajectory a significant portion 
of this change can be provided aerodynami cally (reference 17), with a 
corresponding decrease in fuel requirements. Controlling AMOOS during 
the atmospheric pass of such a mission requires bank angle modulation. 

For instance, the vehicle could fly at the angle of attack for maximum 
lift and a bank angle of 90°. This would give the greatest plr.ne change 
capability and the largest possible lift vector available to use as a 
control. The bank angle could be modulated up or down to correct for un- 
certainties. This would result in a reduced plane change capability 
which could be accounted for in the targeting of the nominal trajectory. 
Since the mission in this study includes a 28.5° plane change, this 
mode was selected for more detailed study. Moreover, this paper will 
develop a guidance law whicu will give only the control history required 
rather than how the vehicle will achieve a particular bank angle. 


ttEPBlODUCKItiW OF 
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SENsiTivrry of amoos to the atmosphere 


The sensitivity of AKOOS to variations in atoos^eric properties 
vas determined hy establishing a ncmtineil trajectozy using a mean atmos- 
phere, generating other trajectories using deviations tram the mean 
atmosphere, and comparing the res\ilts, all other parameters being equal. 

Tvo approaches vere taken to generate the various atmospheres; ccmstant 
multipliers on density and statistical variations of density-altitude 
profiles . 

All of the AMOOS trajectories were calculated with POST, the Pro- 
gram to Optimize Simulated Trajectories (ref. 19). This is a general 
perfomance analysis program which has modularized vehicle and planet 
models and a generalized targeting and optimization capability. It has 
been used extensively on a variety of problems and includes recent 
developments in nximerical integration and optimization techniques. 

The nominal trajectory was calculated using the nominal density 
profile shown in Table 1, the angle of attack corresponding to maximum 
lift, and a bank angle of 90*^ to maximize plane change capability. The 
trajectory program, POST, was used to cecLciilate the magnitude and angle 
of the retro burn at geosynchronous altitude which resulted in the desired 
orbit apogee altitude and inclination upon exit from the atmosphere. 

An iterative minimum norm technique was used for targeting. In this 

case, where the nximber of control veu*iables equals the number of constraints, 

this reduces to the classical Newton-Raphson technique for finding the 
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minisuK of a function. Trajectory Bensitivlties were numerically deter- 
mined \i8ing fcrvard differences based o. ' urbing the control variables. 

To achieve an acceptable phasing orbit, AMOOS must fixit vith an apogee 
between 500 loa and 720 km (reference 17) and an inclination of 28.5^. 

The nominsT was targeted to an apogee of 600 km and an inclination of 28.7^ 
so that ^en a guidance lav S eventually ixsplemented, any remaining slight 
dispersions from the nominal could be tolerated, lu pcurticular, any 
deviations from the nominal bank angle will result in a lovei- value for 
inclination. The values of the trajectory state variables at the point 
of atmosphere entry, 120 km altitude, were used as initial conditions for 
all subsequent trajectories. 

A series of trajectories was run using different ratios of density 
encountered to mean density. This ratio remained constant tfc’Mughout 
each particuleu* trajectory. The result vsis a very wide range of apogee 
altitudes, as shown in figure 7. This figure indicates that very slight 
variations in density (< 13() will result in final apogees outs the 
acceptable range. Density ratios greater than l.Ok result in reentry. 

For densities only slightly lover than the mean, AMOOS exits the atmo- 
sphere vith too much velocity to achieve an acceptable phasing orbit. 

Figure 8 shows the variation of the other critical end condition, orbit 
inclination, vith density ratio. As in the case of final apogee, the 
range of values here is also very vide. Thus, the amount of velocity 
addition, and its corresponding fuel penalty, required to correct for 
errors in inclination is very high. This is shown in figure 9 as a 
plot of velocity addition required to bring the inclination back to the 
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desire-’ value of ?8.5® versus the amount of inclination correction re- 
quiired at various orbit altitudes. 

Constant density ratios throughout a range of altitudes would never 
occur in nature. Therefore, a more realistic atmosphere variation was 
simulated using a statistical model of the atmosphere, developed in un> 
published NASA research, which models over 6,000 Meteorological Rocket 
Network soundings of the atmosphere. The model is in the form of a com- 
puter program which uses a random number generator, assuming temperature 
and density are normally .''.istributed, to generate temperature, density, 
and pressiure profiles whose statistical properties and vex^ical gradients 
match those of the data. The program also provides the mean, standard 
deviation and interlayer correlations of the data. This study used atmo- 
spheres generated for the spring seeison in a latitude band of 15® on each 
side of the equator. Atmospheric data tables were autcanatically input 
to POST from the atmosphere generator program. 

Of the 300 trajectories calculated, each with a different atmosphere, 
137 did not achieve the exit altitude. The final apogees of the remain- 
ing trajectories are presented as a histogram in figure 10. This rein- 
force. the conclusion of the pre-vious figure, that AMOOS is extremely 
sensitive to u.ncertainties in the properties of the atmosphere. 

Further analysis of the state peurameters of the trajectories re- 
vet . .‘d relationships that proved useful in developing a linear feedback 
guidaice lav. Figure 11(a) shows the dependence of the final apogee on 
the vehicle’s orbital energy at exit. Since the vehicle is in an 
elliptic jrbit, the values of energy are negative. The curve is smooth 
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throughout the entire range of trajectories. Figure U(u) is a magnifi- 
cation of the same data in the region of the desired apogee. The allow- 
able apogee range and the corresponding energy range are indicated. The 

2 2 

energy at exit may vary within a band of 300,000 m /s and still result 
in an acceptable trajectory. 

Now, the relationship between final energy and energy during the 

atmospheric pass must be determined. Figure 12 shows the value of energy 

at the point of perigee versus final energy in the region of interest. 

The band of allowable values corresponding to the final energy band 

2 2 

determined in the preceding figure is shown. The 200,000 m /s width 
of the band can be considered acctirate due to the smoothness of the data. 

Similar: results are shown for orbit inclination in figure 13. This 
plot of final inclination versus final energy is also smooth and, along 
with the results of the previous two figures, gives confidence that 
energy management during an AMOOS pass through the atmosphere would lead 
to acciirate targeting of the trajectory. 
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GUIDANCE LAV DEVELOPMENT 


Baised on the previous ANCX)S guidance discussion, a linear feedback 
lav for bank angle modulation was chosen for further study. Two approaches 
were taken to develop this law. First, a relatively simple heuristic 
lav was used to demonstrate the concept. Then an analytical development 
followed using optimal control theory. 

The Eeuristic Approach . The first approach was to have the energy 
of the vehicle follow the energy profile of the nominal trajectory using 
bank angle as the control. The feedback equation was simply, 

non cic b nom 

The nominal bank angle, ^ , was a constant 90** and the gain, K, was a 

noiD 

constant chosen by trial and error with no atteapt at optimization. The 
desired value of energy, E ^^^ , at each point in the trajectory was taken 
from an input table of values obtained from the nominal trajectory. In 
addition, to avoid possible large excursions in bemk angle and resulting 
loss in inclination, the bank angle was limited throughout the treijec- 
torles to a minimum value of 30 and a mAvtnmm value of 130 degrees. 

Various schemes for fumrulating the nomineJ. energy table were tried. 
The best results were obtained by Inputting energy versus flight path 
angle. All trajectories that exit the atmosphere must have a perigee 
point (where flight path angle is zero) and it was shown earlier, in 
figure 12, that matching values of energy at perigee would be advantageous. 
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Also, the initial flight path angle vas Identical for all trajectories 
and the final flight path azkgles for all the trajectories with an accepta>- 
hie final apogee vere nearly the same. Therefore, matching a nominal 
energy versus flight path emgle profile results in an acceptable trajec- 
tory. 

The performance of the heuristic guidance lav is illustrated in figure 
ih. For the range of density factors used the values of final apogee 
generally fell within the acceptable tolerance. The variation of final 
orbit inclination is shown in figure 15. The res\ilt here is excellent 
with the inclination never being more than 0.2 degrees from the desired 
value of 26 . 5 **. Using the results of the last two figures, the level of 
performance of the heuristic linear feedback guidamce lav successfully 
demonstrated the feedback technique and was Judged sufficiently high to 
proceed with an analytical development using optimal control theory. 

The Analyticeil Approach. Since the equations of motion for the atmo- 
spheric portion of AMOOS* flight are highly nonlinear, as is typiccd. with 
entry or skip types of trajectories, the development of an "exact” 
explicit feedback guidance scheme is not feeisible. Therefore a pertur- 
bation guidance scheme was developed involving trajectories in the im- 
mediate neighborhood of the nominal path. By the proi>er choice of qu. \- 
ratic performance criteria and quadratic constraints, a feedback lav will 
be synthesized to keep the system within an acceptable deviation from the 
nominal using acceptable amounts of control. 

Following the analysis of Bryson and Ho (ref. 20), the developnent of 
the optimal regulator is that of the terminal controller with infinite 
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final tine. The eqviatlons of notion of a general system are, 

4 

X«F(X, U), (2) 

vhere X is the n-component state vector and U is the n-component control 
vector. The result of a linear perturbation about a nominal is, 


6X 


AY ^ 


An 


( 3 ) 


For a simpler notation, 

let X = 6X, u = 6U, A « -gy , and B ® ^ 
where A is an n x n and B an n x m sensitivity matrix, 
then (3) becomes, 

X » Ax ♦ , ( 5 ) 

a linear system. To achieve an acceptable level of the perturbed state, 
X, with a minimum of control correction, u, the scalar performance index. 


J = 1/2 ♦ 1/2 j (x^qpt ♦ u^Ru) dt. 


( 6 ) 


o 

is maintained, where S^, Q, and R are positive definite matrices (ref. 20). 
In particular, S^, Q and R are diagonal and, 

s“ = maximum acceptable value of dlag x (t^) x (t^). 


••1 

Q s maximum acceptable value of dlag xx , 

R = maximm acceptable value of dlag uu . 

The scalar Hamiltonian takes the form 


H 





( 7 ) 
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The control, u, that ainiaizes (6) Is obtained bj solving (?) sisoltane- 
ously vith the Euler-’Lagrange equations. 


X • , X(t j) ■ SjX(t j) , 

(8) 

!!■»• 

Performing the differentiations gives. 

(9) 

2. — 
A » -A X -Qx 

(10) 

and -t- Ru ■ 0 ^ u > -R”^B?A . 

(11) 

Substituting (U) into (?) and repeating (10) gives the linear tvo-point 
boundary-value problem, 

X ■ Ax-ffi"S^ x(t^) given (12) 

X » - (^c-A^ X(t^) ■ S^(tj) 

(13) 


In the discussion of the solution of the problea, Bryson and Ho show that, 
for a continuous feedback lav. 



X ■ Sx 

(U») 

Substituting 

ilk) into (12) and (13) gives 



— — -IT — 

X ■ Ax - BR "^B Sx 

(1?) 

and 

(Sx) ■ - A^^. 

( 16 ) 


Performing the differentiation in (l6); 
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REPBODUGISIUTY OF THE 
OBKRNAL PACE IS POOP 



Sx ♦ Sx » -A^Sx - Qx . 


(17) 


Substituting (15) into (17) yields 

Sx + SAx - SB r'^B^Sx * -A*^^ - 5pc, (l8) 

and 

(S + SA - SBR“S^S ♦ A^S + Q) X = 0. (19) 

Since, In general x 0, 

S + SA - SBR’^B^S + A^S + Q = 0. (20) 

Solving for S gives, 

S » -SA - a’^S + SBR"^^S - Q, (21) 

\diich is the oetrix Ricatti Equation. 

Once S is determined, it can be substituted into (lU), vhich in turn 
is substituted into (11) giving 

u « -R“S'^Ebt (22) 

—1 T 

Nov, let C « R B S and substitute into (23) yielding 

u ■ -Cx , (23) 

the desired linear feedback lav. 

A regulator is a feedback controller designed to keep a stationary 
system vithin ea acceptable deviation from a nominal condition. For a 
stationary system, the matrices A and B are constant. Also, let the 
matrices Q and R be constant and consider the case of t^ «. Then the 
matrix Ricatti equation has a possible steady-state finite solution. 


22 



That l8» 

S « 0 » -SA - A^S + SBR'^B^S - Q gives (2l») 

.) -► S at •* «. (25) 

where S is the solution of (2^). 

Then the feedback gain aatrix is 

^ « R"^ ?, (26) 

and 

u » -ac . (27) 

Solution of the Matrix Ricatti Equation * If the system defined by (5) is 
completely controllable and observable, then there is a unique positive 
definite solution. S’, to (2U). (See reference 22.) There are several 
techniques for solving the quadratic equation (2U). For this study, a 
theorea is used to reduce it to a sequence of linear equations, whose 
solutions converge to This technique was presented by deinman in 
reference 23. Equation (2)() is rewritten as 

SA ♦ A*^S - SBR'^'^S ♦ Q » 0. (28) 

Q 

The first step in the procedure is to select a matrix such that A + B T^ 
has all its eigenvalues with negative real parts. Then a sequence of 
matrices S^, S^, 82* • • . , is formed irtiere the S^ are the solutions of 
tue linear algebraic equations 

A^ ♦ a’^ + CjRS^ + Q « 0, (29) 

where 


23 



(30) 


Ci^l - BR 


-1 


and 


- A ♦ B 


' S, « S (reference 2U). 


(31) 


The procedure falls into two phases. First the matrix must he 
found. If A has all its eigenvalues with negative real parts then = 0 
satisfies the requirement. If not, computation of is r< Hi. 

However, since B is a vector in this study, the calculation is 

V 

straightforveurd. In the second phase, equation (20) is solved. Solution 
procedures for such equations are well known and the calculations are 
Sq, ^ quadratic (Reference 23). 

Application to AMOOS Guidance . In POST, the trajectory program, the 
state vector, X, is deflnc»d in a rectanguleu* earth-centered inertial 
reference frame. For convenience, the state vector in this study will 
he defined in the same manner. The state vector has three position com- 
ponents and three velocity components. Since the method is to change 
the cLLtitude profile to control energy, the deviations in the position 
components will not he controlled. The variations of energy seen in the 
section on the sensitivity to the atmosphere result principally from 
changes in velocity. Therefore, the three velocity components of X are 
the only ones of concern to the feedhac.; lav. That is, in this appli- 
cation, the state vector used in the anedytlcal development is assumed 
to contain only the velocity components. Also, as discussed in the sec- 
tion on AMOOS guidance philosophy, the control vector has only one com- 
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poneat, 'bank angle. Therefore, the gain matrix, e , reduces to a vector. 

In the analytlcauL development of the regulator, the A and B matrices 
were assumed constant so that a steady state solution for the gain matrix 
could be found. However, in a typical AMOOS trajectory, the A and B 
matrices change drastically and rapidly during a pass through the atmo- 
sphere. The approach to the problem was then to choose several points 
on the nominal flight path, evaluate A emd B, choose appropriate Q and R 
matrices, solve the matrix Bicatti equation for S, and then calcxilate 
the resxxlting gain vector. The results at all of the points were then 
combined to produce a history of the gain /ector throughout the trajectory. 

Five points along the nominal path were chosen and the rationale for 
their selection is as follows. First, the point of perigee '»as chosen 
since it is that point at which the effects of an unpredictable atmo- 
sphere are greatest and the ability of ■'.he vehicle to control its path 
is the highest. Two points were then selected which essentially marked 
the boundaries of the region in which the vehicle has effective control. 
One of these points was shortly after atmospheric entry, while dynamic 
pressure was increasing, at the threshold of significant aerodynamic 
forces. The other was chosen to occur at a corresponding degree of aero- 
dynamic force, while dynamic pressure was decreasing neeo* the point of 
exit. The final two points were chosen to be halfway between the two 
bovindary points and the point of perigee. Thus, the entire flight spec- 
tznim was included. These five points, the corresponding state variables, 
and other important traject ry parameters are listed li Table 2. 
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The A and B matrices at eech of the five points were calculated 
nunerically, using forward differencing. As a check, scae of the com- 
ponents were re-calculated using central differencing and various per- 
turbation step sizes, with essentially the same results. Each of the 
velocity conponents of the state vector and bank angle were perturbed 
independently and the resixlting changes in F in the questions of notion 
were calculated. The A and B natrices for the z»ainal trajectory are 
presente u. in Table 3. 

The Q and B matric^:> were then calculated as previously described: 

Q - maxinnan acceptable value of diag xx 

-1 — T 

R = BaximuE acceptable v.J.ue of diag uu . 

The process of selecting "acceptable" values of xx and uu was a com- 
bination of two factors. The first involved analyzing the results of 
the atmospheric sensitivity study and the heuristic development. The 
second was to calculate gain vectors using several values for Q and R, 
testing these with a small number of different atmospheres, and ob- 
serving their performance. 

The maxiatum allowable deviation in bank angle finally chosen was 
30°, the same as that used in the heuristic lav. Thus, 

® “ 900 * 

The Q matrix used was detezmined by assuming that the mavlimm deviation 

2 2 

allowed in energy was 100,000 m /s throughout the trajectory. This was 
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translated into allovaUe deviations in the velocity components as 
follows . The equation for energy per unit maas is the sum of the 
kinetic and potential energies and is given hy, 

E = 1/2V^ - ^ , (32) 

r 

where V is the magnitude of the velocity vector, r is the magnitude of 
the position vector, and u is the gravitational constant. This form of 
the energy equation ass\Jones potential energy is zero at infinite radius. 
As discussed previously, the effect of changes in r is neglected and the 
change in energy is given hy. 



Then, 


AE * VAV. 


(33) 

(3l») 


But 


,.2 





1/2 

) 


(35) 


and 

V&V = AV^ + Vg AVg + AV^ , (36) 

where V^, V^, and are the velocity components of the state vector. 

It vas assumed that each component contributed equally to deviations in 
energy. Therefore, at each point in the trajectory, 

V. AV. = •^ or (37) 

1 1 3 

AV. = ^ i = 1, 2, 3 (38) 
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1 


Then, 



(AV^)^ 


« 0 


1 = J 

i,i » 1. 2. 3 

i ^ J 


(39) 


The Q matrix for each of the fire trajectory points is presented in Table 

k. 

Vrosa. reference 20, the systaa x - Ax Bu is completely controllable 

2 

if the rank of the composite matrix (B, AB, A B) is 3. The A and B 
matrices at each of the fire points were used to generate the coaoposite 
matrix euid at eeu:h point the rank of the matrix was three. Thus, the 
systa is con^etely controllable. The system is completely observable 
if there is a matrix D such that 


DD 


Q. (Reference 21 ) 


iko) 


The matrix Q is diagonal and is defined in (39) es* 


1 



» 0 

The matrix D, then, is also diagonal 



i » J 

i,J = 1,2,3. 

i # J 

and can be defined as 

* J 

i,J = 1.2,3. (Ul) 

. # J 


Therefore a matrix which satisfies (UO) exists and the system is completely 
observable. The controllability and observability requiranents are satis- 
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fled and the matrix Ricatti equation has a unique positiv definite 
solution. 

This solution vas found and the gain vector at each point vas cal- 
culated using a RASA developed coopixter program. All of the A matrices 
had eigenvalues vith negative real parts, thus the first phase of the 
Ricatti equation solution vas omitted. All of the solutions converged 
rapidly and the resulting gains are presented in figure l6. Since all 
three curves were smooth, no more trajectory points were required; the 
values of gain between the points could be determined by interpolation. 

The guidance law was synthesized into the trajectory program accord- 
ing to the following block diagram. 



The nominal trajectory and the gain vector were stored versus flight path 
angle as tables vith cubic interpolation, while the nominal control vas 
constant. The bank angle again had fixed limits of 50 and 130 degrees. 
For this einalysis, perfect knowledge of the state vas assioied. 

The feedback guidance lav developed vas tested on a number of atmo- 
spheres provided as before vith the statistical atmosphere generator 
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program. The performance of the lav is illustrated by the histogram in 
figure 17 . Out of the 50 trajectories run, U failed to exit from the 
atmosphere. Kearly all of the final apogees vere lower than desired. 
However, about T2% fell within a 200 Km tolerance, and about 91% vere within 
a 300 Km spread. The mecui of the U6 apogee values s 260.6 and the 
standard deviation is 125.6 !^. I addition, all of the values of final 
inclination for these trajectories were between 27.03'' ur*d 27.91®. -ue 
bank angle histories of these trajectories tended to switch from one fixed 
limit to the other, as illustrated in figure I 6 . The behar'or of this 
particular history is typical. This is most likely due to the difficulty 
of representing the nominal trajectory in a meaningful manner in the 
chosen reference frame. The parameters that can describe a trajectory 
vary widely near the point of perigee throughout the range of atmospheres 
used and, except for energy, there is no partic\ilar reason to constrain 
any of than. Therefore, representing the desired state vector as a func- 
tion of one of the parameters may not have completely satisfactory 
results. 
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CONCLUSIONS MD RECOMffiNDATIONS 


The Aeromaneuvering Orbit-to-Orbit Shuttle (AMOOS) concept for a 
reusable upper stage to be carried by the Space Shuttle has been shown 
to be highly sensitive to \inpredic table variations in the density of the 
atmosphere throug:: which it flies. This was done using a large number 
of atmosphere variations whose statistical properties match those of 
actual data. If the vehicle is not able to compensate for these vari- 
ations, it may fail either to exit from the atmosphere or to dissipate 
enoxigh energy to establish an acceptable phasing orbit with the Space 
Shuttle. Therefore, a guidance law that can compensate for the vari- 
ations is crucial to the feasibility of the AMOOS concept. 

A heuristic linear feedback guidance law* was developed which demon- 
strated that such a law could control the trajectory within acceptable 
tolerances. This law was designed to hold the energy of the vehicle to 
the energy profile of the nominal trajectory. It was tested on a momber 
of different atmosp’^eres using constant multipliers to provide density 
variations. 

An optimal regulator was synthesized using modern control theory. 
The resulting feedback law, based on linearized equations, of motion and 
a study-state assumption, demonstrated the feasibility of this technique. 
This law was tested using the statistically varying atmospheres and the 
values of final orbit apogee were generally lover than the desjred value 
but within a 300 Km band. Also, there was a tendency for the bank angle 
to switch from one fixed limit to the other dwing a trajectory. 
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Several refinements cein be made to the technique to improve perform- 
ance and ease the bank angle switching. One is to define the state vector 
in a more natural reference frame, allowing direct control of important 
parameters such as total velocity. A likely candidate for such a ref- 
erence frame is one which defines the position \'ector by altitude, lati- 
tude and longitude and the velocity vector by velocity magnitude, flight 
path angle, and heading angle. Another advantage cf this system is that 
the sensitivity matrices can be easily determined explicitly rather than 
ntmierically. An even more direct approach would be to use energy itself 
as a state variable (reference 2U). This would have the additional 
benefit of reducing the dimensionality of the system. 

The limitations of linearization could also be reduced. Not only 
is the physical system natvirally non-linear but there is a random term 
introduced by the uncertain variations in the properties of the atmo- 
sphere, A continuous random process term could be added to the equations 
of motion. This term could also be used to account for the absence of 
higher order terns when the system is linearized. The analysis of such 
a system is discussed in tutorial fashion in reference 25. 

The feasibility of AMOOS depends absolutely on its ability to com- 
pensate for unpredictable density variations encountered during the atmo- 
spheric portion of the flight. A feedback guidance lav which provides 
this crucial caj)at*’lity has been developed. Possible refinements to the 
development of the law to further improve its performance have also been 
proposed. 
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TABLE 1. PROPERTIES OF THE MEAN A!IMOSPH£RE 


Altitude, Kin 

Temperattire, ®K 

2 

Pressure, •N/m 

Density 

60 

260 . 70 k 

25.03130 

3.3l*l*8xl0“^ 

63 

255.631 

I6.9621U 

2.3116 

66 

2U1.549 

11.17573 

1.6118 

69 

210. 87U 

6.99571 

1.1557 

72 

208.266 

u. 31921* 

7.221*9x10“^ 

75 

20U.926 

2.61*073 

1*. 1*892 

78 

198.660 

1.69238 

2.9678 

81 

193. 23U 

1.01280 

1.8259 

6 k 

188. UlU 

.59281* 

1.0961 

87 

181. UU2 

.3591*5 

6.9015x10“^ 

90 

177.810 

.20331* 

3.9839 

93 

171.758 

.111*67 

2.3258 

96 

182.265 

.06351* 

1.211*5 

99 

188.930 

.03618 

6.6713x10"*^ 

102 

201.930 

.0217? 

3.7592 

105 

216.930 

.01358 

2.1808 

108 

231.930 

.00871* 

1.3128 

111 

251.930 

.00581 

8.031*1x10“® 

llU 

281.930 

.001*01 

U.9550 

U7 

3U.930 

.00288 

3. 2161* 

120 

3>H.930 

.00213 

2.1701 
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TABLE 2. TRAJECTORY PARAMETERS AT CERTAIN 
POINTS IN THE NOMINAL FLIGHT PATH 


Time, s 

Altitude , 
Km 

Velocity, 

m/s 

Flight Path 
Angle, Deg. 

Energy, 10^ 

2, 2 
m /s 

V^,m/s 

Vgjm/s 

V2,m/8 

0 

120 

10313.8 

-It. 57707 

-8.15lt 

-9532.6 

- 810.5 

-3853.2 

50 

87.371* 

10326.1 

-2.65267 

-8.336 

-91*83.2 

-1281.3 

-3880.2 

90 

73.779 

1005l*.7 

-1.16388 

-11.232 

-9125.8 

-1608.8 

-3902. U 

130 

70.076 

9259.0 

.00553 

-16.951 

-8193.7 

-I81U.O 

-3911.8 

190 

71*. 733 

8337.9 

.87191 

-27.011 

-7076.1 

-2133.7 

-3859.7 

250 

83.810 

8053.0 

I.l8h7 

-29.260 

-6609.1* 

-2589.3 

-3803.1 

»*31 

120 

7950.3 

1.60555 

-20.7Ult 

-5802.2 

-1*098.2 

-3570.3 
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TABLE 3 . THE A ABD B MATRICES AT CERTAIN 
POINTS ON THE FLIGHT PATH 


TIME, S 

* 9 X 

B-^ 
® 30 

50 

-.0003506 

-.C 000291 

.0000909 

-.0000222 

-.0001793 

.0000132 

-.0001592 

-.0000195 

-.0001720 

.0020488 

-.01437 

-.0001894 

90 

-.0029458 

-.0002750 

.0007453 

-.0002530 

-.0015348 

.0001343 

-.0013612 

-.0002318 

-.0014670 

.020855 

-.117939 

.0000731 

130 

-.0047873 

-.0004970 

.0011226 

-.0005259 

-.0025692 

.0002438 

-.0023296 

-.0005065 

-.0024711 

.037871 

-.177564 

.0029141 

190 

-.0019410 

-.0002363 

.0004139 

-.0002941 

-.0011034 

.0001145 

-.0010043 

-.0002894 

-.0010606 


250 

-.0004407 

-.0000622 

.0000918 

-.0000872 

-.0002668 

.0000313 

-.0002313 

-.0000835 

-.0002489 
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TABLE k. THE Q MATRIX AT CERTAIN 
POINTS ON THE PLIGHT PATH 


TIME, 8 

Q MATRIX 

50 

.0809U 0 0 

0 .001U776 0 

0 0 .013550 

90 

.07U952 0 0 

0 . 0023292* 0 

0 0 .013706 

130 

.060U2 0 C 

0 .0029615 0 

0 0 .013772 

I 

190 

i 

.0U506 0 0 

0 .00U0972* 0 

0 0 .0132*08 

o 

u\ 

CM 

.03931 0 0 

0 .00603J* 0 

0 0 .01302 
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1. Delivered to circular Orbit by Space Shuttle 

2. Burn to mission transfer orbit 

3. IVlidcourse correction 

4. Burn to mission orbit 

5. Burn to return transfer orbit 

6. IVlidcourse correction 

7. Burn to Shuttle rendezvous orbit 


Figure L 'Schematic oP Propulsive OOS Mission Profile 
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- Liquid Oxygen Tank 



Figure 2.* Schematic of Propulsive OOS Configuration 
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1. Geosynchronous orbit 

2. Burn to transfer to aeromaneuvering orbit 

3. Miocourse correction 

4. Aeromaneuver to phasing orbit plane and apogee 

5. Phasing orbit apogee , burn to achieve phasing orbit perigee 


Figure 3. - Schematic of AMOOS return trajectory 

ki 



k2 


Figure 4. -AMOOS vehicle configuration 


No Inclination change 



Mission Altitude, km 

Figure 5. - Comparison of AMOOS and propulsive OOS payload capability 



Mach number = 10.27 
Reference Area = 15.69 



Angle of Attack, deg 


Figure 6. - AMOOS hypersonic lift and drag coefficients 
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Figure 7. - Variation of final apogee with density factor 




Inclination correction , deg 

Figure 9. - Velocity addition required to correct inclination error 
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Trajectories 



Apogee at exit , km 


Figure 10. - Histogram of finai apogee 
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Figure 11 (a). - Variation of final apogee with final energy 
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Figure 11 (b). - Magnified view showing acceptable apogee variation 
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Energy at Exit , m^/ $2 


Figure 13. - Variation of final orbit inclination with final energy 
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Figure 13. - Variation of final inclination with density 
factor for the heuristic guidance law 
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Figure 17. - Histogram of final apogee for optimal regulator test 
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Figure 18. - Typical bank angle history from optimal regulator test 
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